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Abstract 

The 19 F(a,p) 22 Ne reaction is considered to be one of the main sources of fluorine depletion in 
AGB and Wolf-Rayet stars. The reaction rate still retains large uncertainties due to the lack of 
experimental studies available. In this work the yields for both exit channels to the ground state and 
first excited state of 22 Ne have been measured and several previously unobserved resonances have 
been found in the energy range E/ a 6= 792- 1993 keV. The level parameters have been determined 
through a detailed R-matrix analysis of the reaction data and a new reaction rate is provided on 
the basis of the available experimental information. 

PACS numbers: Valid PACS appear here 
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FIG. 1: Energy level scheme for the 19 F(a,p) 22 Ne reaction. The entrance channel (a+ 19 F) has a 
threshold of 10.47 MeV with respect to the ground state of 23 Na. The compound state can decay 
either via the po channel to the ground state of 22 Ne or by the pi channel to the first excited 
state of 22 Ne. The subsequent emission of a 1.27 MeV photon from the decay of the first excited 
state (2+) to the ground state (0+) of 22 Ne can also be observed. The dark rectangle above 23 Na 
represents the energy region studied in this work. White regions above and below it correspond to 
energies where the cross sections were extrapolated. 

I. INTRODUCTION 

Fluorine is by far the least abundant of the elements with atomic mass between 11 and 32. 
While several nucleosynthesis scenarios have been proposed for the production of fluorine, a 
unique site for the origin of this element has not been identified yet. 

Presently three different scenarios are being discussed for the origin of fluorine. One 



Q. 



A second 



possible process is the neutrino dissociation of 20 Ne in supernovae type II 
scenario is the pulsating He-burning stage in AGB stars |2| and the third possibility is the 
hydrostatic burning of helium in Wolf-Rayet stars [3j . It may be possible that all three sites 



contribute to the formation of fluorine in the universe 



Q. 



So far, the on 



objects in the galaxy where fluorine has been observed are AGB stars [5] and post-AGB star 
configurations |g|. 

Fluorine nucleosynthesis in AGB stars takes place in the hydrogen-helium intershell region 
where the 14 N ashes from the preceding CNO burning are converted to 18 F by a-particle 
captures (a similar reaction path is followed by Wolf-Rayet stars). The unstable 18 F isotope 



y extra solar system 
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decays with a half life of 109.8 minutes to 18 0. Both a proton or an a-particle can be 
captured by 18 0. In the former case, 4 He and 15 N are being produced, while in the latter 
case 22 Ne is formed. This second possibility does not produce fluorine. The alternative 
18 0(p,a) 15 N(a, 7) 19 F capture reaction is the main fluorine production link in this scenario. 
The abundance of 15 N depends sensitively on the hydrogen abundance in the inter-shell 
region. Hydrogen has been depleted in the preceding CNO burning but can be regenerated 
through the 14 N(n,p) 14 C reaction, with the neutrons being produced by the 13 C(a,n) 16 
reaction. Additional protons may be mixed into the region when the convective envelope 
penetrates the intershell region at the end of the third dredge up (TDU). 14 C produced by 
this reaction provides a second link for the production of 18 via 14 C(a,7) 18 0. 

Fluorine is very fragile and three reactions may cause rapid fluorine destruction. Because 
of the high abundance of 4 He in the intershell region, the 19 F(a,p) 22 Ne is expected to be a 
dominant depletion link. Another depletion process corresponds to the 19 F(n,7) 20 F reaction, 
caused by neutrons being produced by the 13 C(a,n) 16 or the 22 Ne(o;,n) 25 Mg reactions. If 
hydrogen is available in sufficient abundance, fluorine may also be depleted through the 
very strong 19 F(p,a) 16 reaction. All this, however, is strongly correlated with the existing 
proton and a-particle abundance at the fluorine synthesis site. 

The reaction rate for 19 F(a,p) 22 Ne is highly uncertain at helium burning temperatures. 
Even recent nucleosynthesis simulations ?! still rely on the very simplified rate expression 
of Caughlan and Fowler (CF88) [8] based on an optical model approximation for estimating 
the cross section of compound nuclear reactions with overlapping resonances. This ap- 
proach was originally developed [s| and employed 1(3, H] for cases where no experimental 
information was available. This reaction rate is in reasonable agreement with more recent 
Hauser-Feshbach estimates assumin g a high density of unbound states in 23 Na [12]. Other 
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19 F nucleosynthesis simulations 

19 F(a,p) 22 Ne for Ez a 6=1.3 MeV to 3.0 MeV [151 ]. estimating also possible low energy resonant 



14] rely on the sparse experimental data available for 

2 possible low energy resonant 
16]. For example, in their ap- 



contributions from known a-particle unbound states in 23 Na 
proach they approximated the reaction rate by only considering the resonance contributions 
of the various states neglecting possible interference effects and broad-resonance tail contri- 
butions. In this paper we describe a new measurement of 19 F(a,p) 22 Ne at lower energies. 
In our study we explored the reaction for E; a &=792 to 1993 keV at different angles. The 
pronounced resonance structure was analyzed in terms of the multi-channel R-matrix model 
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using the recently developed R- matrix code AZURE [17] . Energy regions not measured here 
were also included in the computation of the stellar reaction rate by using a combination 
of data available in the literature and our experimental results to extrapolate the reaction 
cross section. Possible low energy resonances were taken into account sampling characteris- 
tic a-particle partial widths of the known a-particle unbound states in 23 Na, while deriving 
the corresponding proton partial width from the available elastic scattering data in the lit- 



erature 



18]. Interference terms were modelled with Monte Carlo simulations. In section 2 



we describe the experimental setup and the experimental procedure. The R-matrix analysis 
of the experimental data and the determination of the nuclear parameters of resonances is 
described in section 3. Finally, the experimental and extrapolated reaction rates (with error 
bars) are presented and discussed in section 4. 



II. EXPERIMENTAL METHOD 

The destruction of 19 F by an a-particle capture occurs mainly by a resonant reaction 
process through the 23 Na compound nucleus in an excitation range of high level density. 
The populated resonant levels will decay by proton emission to the ground state (po) or first 
excited state (pi) of 22 Ne. This study included the direct measurement of the 19 F(a,p) 22 Ne 
reaction by detecting both p and pi protons (see 19] and [^(J) as well as the additional 
measurement of the pi channel via the detection of the 7-ray transition 2jJ from the first 
excited state (2 + ) to the ground state (0 + ) of 22 Ne (see figured]). 

The experiment was performed at the Nuclear Science Laboratory of the University of 
Notre Dame using the 3.5 MV KN Van de Graaff accelerator. In a first run, the excitation 
curve of 19 F(a,p) 22 Ne was investigated between Ei a b=122A keV and 1993 keV. For this 
experiment two Si surface barrier detectors with a depletion depth of 300 /^m were mounted 
at forward angles, while one 100 /im Si detector was positioned at a backward angle. These 
thicknesses were sufficient for stopping the reaction protons. The effective solid angle of the 
detectors at each position configuration was determined using a mixed 241 Am and 148 Gd 
a-particle source with a known activity placed at the target position. 

The energy range between E; a {,=1629 keV to 1993 keV was mapped in 5 keV energy 
steps with the detectors mounted at 30°, 90° and 130°. This made it possible to use the 
two known resonances at E; a £,=1.67 MeV and 1.89 MeV [ijj] as reference for calibrating and 
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matching the reaction yield to the previous results. At lower energies, the detector position 
was changed to 40°, 100° and 120° with respect to the beam direction. The excitation curve 
was mapped from Ez a 6=1220 to 1679 keV using the E; a ft=1.67 MeV resonance as a reference. 
In total 483 proton spectra were acquired and for every energy, one elastic scattering a- 
particle spectrum was taken. The stoichiometry of targets was constantly monitored with 
back-scattered a-particles. 

The 19 F transmission targets were prepared by evaporating 10 /ig/cm 2 of CaF2 on 40 
/zg/cm 2 natural carbon substrates, mounted on aluminum frames. The target was placed 
with the evaporated material facing the beam on a ladder attached to a rotating rod at the 
center of the scattering chamber. The ladder held one target and a collimator that was used 
for centering the beam. The targets deteriorated significantly under beam bombardment, so 
their stability was monitored by measuring frequently the yield of the elastically scattered 
a-particles at the E iafe =1.89 MeV resonance. Targets were constantly replaced with new and 
recently evaporated targets. 

For any two detectors with the same absolute efficiency the relative count rate is inde- 
pendent of target stoichiometry and beam intensity, as expressed by 



f da(E,9) y 1 f da(E,e) \ 

iab\ dn ) 2 \ dn )i u 

cm 



Ni(E, 9) and N 2 (E, 9) corresponds to the number of events in detectors 1 and 2, respectively, 
[*§p] iab/cm is the solid angle correcting for center of mass to the laboratory system, and da( ^ e ^ 
are the differential cross sections measured at detectors 1 and 2, respectively. The differential 
cross section of the 19 F(a,p) 22 Ne reaction was determined relative to the differential cross 
section for elastic scatte ring measured at 160°. It has been shown by Huang-sheng et al. 



22j and by Cseh et al. 23j that below E; a ;,=2.5 MeV the elastic scattering of 4 He on 19 F 



follows the Rutherford law, 

da(E,0) da(E,0) ( Z x Z 2 e 2 ^ 2 



, sin" 4 -. (2) 

dVL el dVL Ruth V AE 

Z% and Z 2 are the atomic numbers of projectile and target, respectively, e is the proton 
electric charge and e 2 =1440 keV fm, E is the relative energy of target and projectile in the 
center of mass system, and 9 is the center of mass angle at which the elastically scattered 
particles are observed. Within the small thickness of the target (27±5 keV) the stopping 
cross section e is assumed to be a constant. The variation of the elastic differential cross 



6 



section across the target thickness is also very small. Therefore the elastic yield can be 
expressed as 

Y^(*<m * (3, 

V dVl J Ruth e 

Subsequently the target integrated proton yield Y p can be expressed relative to the elastic 
scattering cross sections as 



Y n =- 



The energy dependence of the stopping power dE / dx of 4 He is well known for both calcium 



A/ da(E,6) \ ( N prot (E,9) \( dQ prot \ A f Z 1 Z 2 e 2 V . _,9 
*V d ^ J, rot \Neias(E,6)){dn ela J lab/cm e{ AE J Sm 2- U 



and fluorine 



241 ]. The stopping power for fluorine was fitted to a quadratic polynomial 
function given by 



dx j F ( i 

A similar relation was determined for the calcium stopping power. The partial stopping 
cross section e for each of the nuclear species in the target is described by: 

1 f dE \ 

where 

n = u P N A /A (7) 

with v the number of atoms per molecule, p the density of the target (again assuming the 
evaporated material has the same density as the powder used before target preparation), 
Na the Avogadro number, and A the mass number. The total stopping cross section of the 
calcium fluoride target depends critically on the target stoichiometry 

i N Ca , v 

e = e F + —e Ca . (8) 

The ratio measured in the evaporated target layer does not necessarily reflect the sto- 



ichiometry of the material before being evaporated. It was reported in previous work [25] 
that evaporated CaF2 shows a stoichiometric calcium to fluorine ratio of 1:1. The targets 
used in the present experiment were tested using well known resonances in the 19 F(p,«7) 16 



reaction. The results indicated 



as that of the CaF 2 powder [26]. 



:hat the stoichiometry of the evaporated material is the same 
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FIG. 2: Proton spectra for 19 F(a,p) 22 Ne at 135° for three beam energies (label at the right upper 
corner of each box). The upper spectrum shows both proton groups at the reference resonance. 
Both peaks appear clean of background and are easy to identify. The spectrum at E; a fe=1100 
keV shows two groups of protons, still well isolated from the background. Finally, a spectrum at 
E; a ;,=792 keV, where no proton groups were positively identified is shown. The integrated charges 
are (in pC) 1981, 180647, and 519252, respectively. 



In a second set of experiments we measured the proton yield at lower energies. The 
experimental setup was designed to perform the measurements with higher beam currents 
and larger detector solid angles to compensate for the drop in reaction cross section. 

The target chamber for this set of experiments allowed mounting of the target at two 
different angles with respect to the beam: at 45° and 90°. With the first option we measured 
scattering angles below 90°, while with the second other angles were measured at a smaller 
effective target thickness. We tested several 19 F-implanted targets for stability. Substrates 
tested were Ta, Ni, Cr, Al, Fe, and Mo. The best stability against beam deterioration was 
obtained for the Fe substrate. Electron suppression was supplied through an aluminum 
plate at -400 volts, mounted 5 mm in front of the target. Carbon buildup on the target was 
minimized with a liquid nitrogen-cooled copper plate. The target itself was water cooled 
from the back and electrically isolated from the scattering chamber. Beam current was 
directly measured at the target holder. 

Two Si detectors were mounted on the rotating plate with aluminum holders. Collimators 
were placed in front of the detectors and pin hole- free Al foils were used to stop the elastically 
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scattered a-particles, while allowing the protons to reach the surface of the detectors. Both 
detector holders and the rotating plate were electrically isolated from the rest of the chamber. 
The detectors had an effective detection area of 450 mm 2 . The solid angles for both detectors 
were determined using a mixed 241 Am+ 148 Gd a-particle source of known activity mounted 
at the target position. We measured them to be 0.130±0.026 and 0.133±0.027 steradians, 
respectively. 

With the chamber at the perpendicular target position and both detectors at 135°, a to- 
tal of 540 spectra were acquired for laboratory energies between 792 and 1380 keV. Typical 
spectra are shown in figure [2} the last spectrum taken (E; a b=792 keV) does not show iden- 
tifiable proton groups. Subsequently, the chamber was reoriented to the 45° target position. 
The detectors were mounted at 150° and 120° with respect to the beam direction and 178 
spectra were taken. Finally, the detectors were mounted at 75° and 105° with respect to the 
beam direction and 69 spectra were measured. 

In this experimental configuration the reaction yield could not be normalized to elastic 
scattering because of the thick target backing material. The yield was therefore measured 
relative to the accumulated charge on the target during each run. The reaction yield Y p (9) 
for the detector mounted at angle 9 is derived from the number of registered proton events 
N p (#) in the detector is 

Y (9)= N ^ (9) 

where N Q is the accumulated charge, e p is the absolute detection efficiency (assumed to be 
1 for charged particles and Si detectors at very low count rates), and dfl p is the solid angle 
subtended by the detector. 

The differential cross section can be directly derived from the reaction yield normalized 
to the yield of the Ej a &=1.37 MeV resonance as measured in the first experiment. This 
depends critically on the stability of the fluorine content in the target material. Since 
the amount of fluorine decreases with accumulated charge, the reaction yield needs to be 
corrected accordingly. During the experiment the yield of the Ej a 6=1.37 MeV resonance was 
monitored frequently to correct for target degradation. The on-resonance thick target yield 
Y tt as a function of accumulated charge N Q can be expressed by the linear relation: 

Y tt = a + bN a , (10) 
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with a and b as constants. The measured yield of the observed protons Y P (E) was corrected 
for target degradation in terms of the accumulated charge to 

Data obtained from all the described experiments consist of 20 excitation functions, with 
eleven corresponding to 19 F(a,po) 22 Ne and nine to 19 F(a,pi) 22 Ne. Ten angles were measured 
in different energy regions. All add up to 1471 data points (See tabled]) which were analyzed 
in terms of the R-matrix theory. 

III. MULTICHANNEL R-MATRIX ANALYSIS 

ysis of the experimental data we used the A-matrix version of the computer 



For the ana 

code AZURE 17|. AZURE is a multichannel and multilevel code that implements the A- 
and R-matrix formalisms as presented by Lane and Thomas {27]. The code is capable of 
fitting experimental datasets by varying formal parameters (energies and width amplitudes) 
of compound-nucleus states. The integrated cross section can also be computed from angular 
distribution datasets. Error bars for both the parameters and cross section curves are treated 
with Monte Carlo techniques. 

The input to the code consists of a set of initial nuclear parameters; each level is char- 
acterized by one energy eigenvalue and several formal reduced width amplitudes (one per 
channel per level). The code identifies open reaction channels and from the input it assigns 
an independent width amplitude for each (s,l) combination allowed for the level. Theoret- 
ical differential cross section curves at different angles are computed and then compared 
to experimental yields after target integration corrections 281 ] . The maximum likelihood 
estimator is then minimized by varying all the parameters simultaneously. Each time a local 
minimum is found, the integrated value of the cross section is computed. 

Overlapping of resonances complicated the simultaneous fitting of the yield curves. The 
code was by itself not able to find a set of formal parameters that would reasonably describe 
the yield curves. For this reason, an initial parameter set had to be found without the help 
of the minimization routines. By trial and error, the choice of initial nuclear parameters 
was done by adjusting the energy for each of the levels so as to describe the position of 
resonances as close as possible. The height and width of the resonances was, on the other 
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TABLE I: Yield curves measured in this work 



curve 


channel 


angle 


Emiv, fkeV) 


YjmnT (ke\0 

77 LUX \ ^ / 


A (keV) 


1 


19 F(a, Po ) 22 Ne 


130 


1641 


1993 


15 


2 


19 F(a,p ) 22 Ne 


90 


1641 


1993 


15 


3 


19 F(a,p ) 22 Ne 


30 


1641 


1993 


15 


4 


19 F(a,p ) 22 Ne 


120 


1224 


1679 


15 


5 


19 F(a,p ) 22 Ne 


100 


1224 


1679 


15 


6 


l9 F(a,p ) 22 Ne 


40 


1224 


1679 


15 


7 


l9 F(a,p ) 22 Ne 


105 


1027 


1367 


25 


8 


19 F(a,p ) 22 Ne 


120 


929 


1359 


35 


9 


19 F(a,p ) 22 Ne 


135 


792 


1363 


25 


10 


19 F(a,p ) 22 Ne 


150 


929 


1359 


35 


11 


19 F(a,p ) 22 Ne 


75 


1027 


1367 


35 


12 


19 F(a vi) 22 Ne 


130 


1629 


1981 


15 


13 


l9 F(a,pi) 22 Ne 


90 


1629 


1981 


15 


14 


l9 F(a, Pl ) 22 Ne 


30 


1629 


1981 


15 


15 


l9 F(a, Pl ) 22 Ne 


120 


1224 


1679 


15 


16 


l9 F(a, Pl ) 22 Ne 


100 


1224 


1679 


15 


17 


19 F(a, Pl ) 22 Ne 


40 


1224 


1679 


15 


18 


19 F(a, Pl ) 22 Ne 


120 


929 


1359 


35 


19 


l9 F(a, Pl ) 22 Ne 


135 


792 


1363 


25 


20 


l9 F(a, Pl ) 22 Ne 


150 


929 


1359 


35 



hand, approximated by varying the width amplitudes. 

Interference patterns between the various resonances were determined by iteratively prob- 
ing the contribution to yield curves of levels within groups of the same P". The sign of the 
width amplitudes was flipped one by one for all the (s,l) channels. These steps were repeated 
iteratively several times until the theoretical curves resembled the dataset. The resulting set 
of parameters was then used as input to the R-matrix code coupled to the x 2 minimization 
routines. Every time a calculation was performed all 20 excitation curves were examined. 
The target thickness A used for each of the yield curves is shown in table [B 
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The resulting R-matrix fits to the experimental yield data Y p are shown in figures [3] and 
H] for the ground state and the excited state transitions, respectively. Both the energies 
and reduced width amplitudes 7, were determined with a single channel radius (a c =5.5 fm). 
The boundary condition was set to S c , the shift function, at the level in each J 71 " group with 
the lowest energy. Background states were included for each of the J 71 " groups. (The set of 
201 R-matrix parameters and the experimental dataset can be obtained by contacting the 
author.) 

Parameter error bars were determined with a bootstrap method and correspond to 
a confidence interval of 95% (see figure [5]). From the set of 1,471 experimental data points 
we generated 40,000 subsets, with each subset containing 1,471 data points as well. The 
selection of points was performed with a Monte Carlo method and as a result, each set 
contains a random number of points repeated more than once. Using the set of formal 
parameters obtained with the R-matrix analysis described above, a x 2 was computed for 
each of the subsets sampled by bootstrapping. The distribution of \ 2 values is shown in 
figure [3(a). The 0.95 cumulative value of the distribution corresponds to x 2 = Xgs = 22,774. 
Finally, with a Monte Carlo technique, we generated x 2 curves by varing each parameter 
independently while fixing all others. The error bar corresponds to the largest parameter 
value such that x 2 < X95 ( see figure [5(b)). Total cross sections for both channels were 
computed as well and are discussed in the next section and shown in figure [61 



IV. THE REACTION RATE 



The thermonuclear reaction rate has contributions from three energy regimes: a) the 
region measured experimentally in this work, spanning from E/ a b=792 keV to 1993 keV, b) 
the region below E; afe =792 keV (not measured), and c) the region above E iafe =1993 keV (not 
measured here as well). 

The contribution to the rate from our experimental dataset was calculated by integrating 
numerically the total cross section a(E) over the Maxwell-Boltzmann distribution of stellar 
gas at temperature T 

N A (av) = Na(JA ' (kT)- 3 / 2 J a(E)Eexp(^jdE. (12) 
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FIG. 3: Experimental yield curves for the 19 F(a,po) 22 Ne channel and the R-matrix fits (solid lines). 

Different curves are offset by a factor of 100. Vertical axis units are relative yields. The labels to 

the right of each curve are described in table [H 
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FIG. 4: Experimental yield curves for the 19 F(a,pi) 22 Ne channel and the R-matrix fits (solid lines). 
Different curves are offset by a factor of 100. Vertical axis units are relative yields. The labels to 



the right of each curve are described in table [H 
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FIG. 5: Calculation of error bars, (a) Bootstrap histogram used to compute the error bars with 
a 95% confidence level. The horizontal axis corresponds to the raw x 2 value associated to the 
set of central-value parameters. Each bin contains intervals Ax 2 =20. A total of 40,000 subsets 
were sampled from the experimental set with a Monte Carlo method. Note that the curve does 
not correspond to a normal distribution. However, the centroid corresponds to the raw value of 
X 2 associated to the best R-matrix fit (Xmm = 20>535). An area corresponding to 95% of the total 
integral is below Xg5=22,774 (vertical line), (b) Sample x 2 vs - parameter value (7) curve. The 
horizontal line represents ^95, while the dotted vertical lines are, from left to right, the lower limit, 
central value, and upper limit of the parameter, respectively. An equivalent curve was generated 
for each of the formal parameters. 
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FIG. 6: R-matrix calculated total cross sections for 19 F(a, p) 22 Ne. Both 19 F(a, po) 22 Ne (upper 
panel) and 19 F(a, pi) 22 Ne (lower panel) curves are shown with uncertainty bands, as computed by 
sampling the 95% confidence interval nuclear parameter space with a Monte Carlo method. The 
strong oscillations in the upper and lower limits of the bands are an artifact of the finite number 
of samples (10,000) taken with the Monte Carlo method. 

Here E is the energy of the particles in the center of mass system, is Avogadro's number, 
k is Boltzmann's constant, fi the reduced mass, and T the temperature of the gas. The 
total cross section a was derived from the R-matrix calculation with the recommended 
values of formal parameters. Both po and pi components are shown in figure EJ Upper 
and lower limits of the total cross section and the reaction rate in this energy regime were 
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computed by sampling the parameter space defined by the upper and lower values of the 
fitting parameters with a Monte Carlo technique. All parameters were varied simultaneously 
and a total of 10,000 parameter sets were produced. The reaction rate for each parameter 
set was calculated with equation [T2] and all resulting rates were compared with each other. 
The highest (lowest) value obtained corresponds to the upper (lower) limit of the reaction 
rate. 

The contribution to the reaction rate from resonances below Ej a b=792 keV, which were 
not measured in the present experiment, was considered as well. Previous studies through 
other reaction channels do indicate several unbound states in 23 Na in this energy range 
near the a-particle threshold [16], which may contribute significantly to the 19 F(a,p) 22 Ne 



reaction. Most notably, detai 
this energy range in 23 Na 



ed elastic proton scattering measurements were performed for 
18] , and provide important information necessary for estimating 
the contributions of these lower energy states to the reaction rate. Resonances observed 
in the 22 Ne(p,p) 22 Ne and 22 Ne(p,p') 22 Ne channels were used to define the spins, parities, 
energies, and both proton po and pi partial widths of contributing states. On the other 
hand, a-particle partial widths r a were obtained by adopting the experimentally known 
a-particle reduced widths 7 2 determined in the high energy range. 

Energies and reduced widths obtained with the R-matrix analysis were used to calculate 
the logaritmic average < log(j^) > for each J 71 " group (see figure [7j). The extrapolated 
reduced a-particle width amplitude for a state with parameters (J,7r) was chosen to be 

7 ^(J,7r) = 10 </o ^ )> . (13) 

Extrapolated upper (lower) limit values of the reduced a-particle widths were set equal 
to the highest (lowest) 7 2 value determined from the experimental data of the corresponding 
J 71 " group. The set of extrapolated 7 2 values, together with the reduced proton widths 7 2 
calculated from 22 Ne(p,p) 22 Ne and 22 Ne(p,p') 22 Ne experiments by r p = 2P7 2 (where P is 
the penetrability through the Coulomb barrier) is shown in table [Til 

Based on these parameters, the resonant contribution to the cross section was calculated 
using the R-matrix formalism. Nevertheless, the interference pattern between resonances 
can not be predicted in this scheme. Therefore, we simulated its effect with a Monte Carlo 
sampling of possible width amplitude sign combinations. Figure M shows the resulting S- 
factor for three sample assumptions for interference between resonances. Upper and lower 
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TABLE II: Nuclear parameters for states with E lab <792 keV. 



J 


7T 


£ cm (MeF) 


ll (MeV) 


ll (MeV) 


ll {MeV) ' 








Recomm 


Upper 


Lower 


1.5 


1 


0.010 


3.72E-4 


5.29E-4 


2.61E-4 


1.5 


-1 


0.031 


8.89E-5 


1.51E-4 


5.37E-5 


0.5 


1 


0.037 


9.60E-6 


2.44E-4 


1.48E-6 


2.5 


1 


0.049 
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FIG. 7: The experimentally determined 7^ as a function of center-of-mass energy, presented in sets 
of J\ 

values of the reaction rate were calculated with equation [12] from different combinations of 
signs and widths sampled within the parameter space, as before. The recommended value 
corresponds to the logaritmic average of the upper and lower limits of the reaction rate. 

Finally, the contribution to the reaction rate from resonances above E; a b=1993 keV was 
computed by extrapolating our experimental rate to higher temperatures by following the 
energy dependence of the Hauser-Feshbach rate MOST 3o|. We did not perform an R- 
matrix analysis due to the current uncertainty of spins and parities of excited states in 
23 Na in this energy region, for which several experimental works have been published (for 



example see 



15 



23 



31 



321 ] ). Over one hundred resonances have been identified for 2.0 < 



Eiat (MeV) < 4.7, but the data has not been able to constrain the spins and parities of most 
of the states. The average level density is 0.04 keV -1 , high enough to apply the Hauser- 
Feshbach formalism for calculating the reaction rate 33J . At the upper limit of the Gamow 
window corresponding to our experimental data (T=lxl0 9 K), the agreement between the 
Hauser-Feshbach and our R-matrix calculated rate is very good. We used the experimentally 
determined reaction rate here to renormalize the MOST values above this temperature. 

The total reaction rate consists, for temperatures below T=lxl0 9 K, of the sum of the 
rate from the R-matrix analysis of our experimental data and the rate calculated from 
extrapolated Monte Carlo cross sections at the lowest energies, and of the Hauser-Feshbach 
renormalized rate above T=lxl0 9 K. The resulting reaction rate is shown in table UTT1 
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TABLE III: Reaction rate for 19 F(a,p) 22 Ne 



n 


Na < crv > recomm 
[cm 3 / 's mol] 




Na < av > low 
[cm? / s mol] 


Na < (tv > up 
[cm? / s mol] 


0.10 


2.402E-22 




1.049E-23 


5.500E-21 


0.11 


5.072E-21 




4.173E-22 


6.166E-20 


0.12 


6.322E-20 




8.649E-21 


4.621E-19 


0.13 


5.625E-19 




1.158E-19 


2.732E-18 


0.14 


3.943E-18 




1.117E-18 


1.392E-17 


0.15 


2.477E-17 




8.120E-18 


7.555E-17 


0.16 


1.399E-16 




4.807E-17 


4.070E-16 


0.18 


2.758E-15 




1.073E-15 


7.091E-15 


0.20 


3.310E-14 




1.455E-14 


7.556E-14 


0.25 


3.680E-12 




2.047E-12 


7.088E-12 


0.30 


1.272E-10 




7.709E-11 


2.759E-10 


0.35 


2.431E-09 




1.484E-09 


6.038E-09 


0.40 


3.340E-08 




1.939E-08 


7.961E-08 


0.45 


3.631E-07 




2.064E-07 


7.438E-07 


0.50 


3.072E-06 




1.796E-06 


5.475E-06 


0.60 


1.020E-04 




6.176E-05 


1.601E-04 


0.70 


1.445E-03 




8.658E-04 


2.151E-03 


0.80 


1.115E-02 




6.630E-03 


1.624E-02 


0.90 


5.615E-02 




3.346E-02 


8.186E-02 


1.00 


4.173E-01 




2.483E-01 


6.068E-01 


1.25 


5.748E+00 




3.398E+00 


8.746E+00 
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FIG. 8: Monte Carlo simulated low energy S-factor for 19 F(ct, p) 22 Ne for three different sample 
resonance interference assumptions. 

Figure [9] shows the total reaction rate relative to the phenomenological rate estimate of 
CF88. Shown are the Hauser-Feshbach model predictions using the codes MOST and NON- 
SMOKER relative to the CF88 predictions. Also compared is the rate estimate of Lugaro 
et al. [ljj, which is based on a single non- interfering resonance approximation. 

The new rate is significantly higher (about one order of magnitude in the stellar tem- 
perature regime) than the rate based on the assumption of single non-interfering resonance 



levels of Lugaro et al. 



131 ] . This can be attributed to the fact that the new rate is calculated 



from non-narrow resonance contributions, as given by the R-matrix analysis. Also, the CF88 
rate is very similar to the new rate, except for the astrophysically relevant temperature T 
= 3xl0 8 K, where it is close to a factor of 10 smaller. However, the statistical model pre- 
dictions overestimate the reaction rate and can be rejected at a 95% confidence level for 2.5 
< T (10 8 K) < 5.5. For the highest temperatures, the new recommended rate differs from 
the statistical model predictions by a renormalization factor (0.62). 

The main source of uncertainty in the new rate at AGB star temperatures comes from 
the uncertainty in the extrapolation of the reaction cross section. Experimental work below 
E/a6=800 keV is required to constrain the partial a-particle widths and the interference pat- 
terns in the 19 F(a,p) 22 Ne reaction. We have shown the importance of the interference effects 
between resonances. Therefore, a direct measurement towards lower energies is probably the 
only plausible solution to the problem of the uncertainty of this reaction rate at AGB star 



21 



100 



~i 1 1 1 — i — i — r 



t 1 1 1 — i — i — r 



OO 

oo 

u 

o 



10^ 




0.1 



> 



0.0 



6 



j i i i i i i_ 



Lower limit 
Upper limit 

Recommended 
MOST 
Lugaro et al. 

NON-SMOKER 

i i i i i i i_ 



1 



10 



T(10 K) 
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131 ] rate. 



le F(a,p) Ne reaction compared to the phenomeno- 
, the Hauser-Feshbach predictions, and Lugaro et aZ.'s 



temperatures. 

For explosive stellar scenarios (T > lxlO 9 K ) the situation is still more delicate as 
spins and parities of resonances contributing to the rate are uncertain. Both direct and 
indirect measurements of the 19 F(a,p) 22 Ne reaction above E/ a (,=2 MeV can help improving 
the quality of the rate for explosive scenarios. 
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V. SUMMARY AND CONCLUSIONS 

We have measured the 19 F(a,p) 22 Ne reaction in the energy range £^6=792-1993 keV. 
Stable fluorine targets were developed and several resonances were found in the 20 experi- 
mental yield curves. Ten different angles ranging from 30° to 150° were measured and two 
reaction channels (po and pi) observed. An R- matrix analysis of the data was performed 
to determine the differential and total reaction cross sections in the investigated excitation 
energy range. The cross section is characterized by many broad resonances tailing into the 
low energy range. Possible additional resonance contributions in that excitation range were 
predicted in a Monte Carlo cross section analysis on the basis of available data on the nuclear 
level structure of the 23 Na compound nucleus. The predicted contributions were included in 
the final reaction rate analysis. A full analysis of the impact of this new rate on the fluorine 
production in AGB stars will be presented in a subsequent paper. 



[1] S. E. Woosley and W. C. Haxton, Nature 334, 45 (1988). 

[2] M. M. Busso, in Planetary Nebulae in our Galaxy and Beyond, edited by M. J. Barlow and 

R. H. Mendez (2006), vol. 234 of IAU Symposium, pp. 91-98. 
[3] S. Goriely, A. Jorrisen, and M. Arnould, On the mechanisms of 19 F production (Max-Planck- 

Institut fiir Physik und Astrophysik Rep., 1989). 
[4] A. Renda, Y. Fenner, B. K. Gibson, A. I. Karakas, J. C. Lattanzio, S. Campbell, A. Chieffi, 

K. Cunha, and V. V. Smith, Mon. Not. R. Astron. Soc. 354, 575 (2004). 
[5] A. Jorrisen, V. V. Smith, and D. L. Lambert, Astron. Astrophys. 261, 164 (1992). 
[6] K. Werner, T. Rauch, and J. W. Kruk, Astron. Astrophys. 433, 641 (2005). 
[7] A. Palacios, M. Arnould, and G. Meynet, arXiv:astro-ph/0508031 (2005). 
[8] G. R. Caughlan and W. A. Fowler, Atom. Data Nucl. Data 40, 283 (1988). 
[9] W. A. Fowler and F. Hoyle, Astrophys. J. Suppl. S. 9, 201 (1964). 
[10] R. V. Wagoner, Astrophys. J. Suppl. S. 18, 247 (1969). 

[11] W. A. Fowler, G. R. Caughlan, and B. A. Zimmerman, Annu. Rev. Astron. Astr. 13, 69 
(1975). 

[12] F.-K. Thielemann, M. Arnould, and J. W. Truran, MPA Rep., No. 262, 15 pp. 262 (1986). 



23 



[13] M. Lugaro, C. Ugalde, A. I. Karakas, M. Wiescher, J. Gorres, J. C. Lattanzio, and R. C. 

Cannon, Astrophys. J. 615, 934 (2004). 
[14] R. Stancliffe, M. Lugaro, C. Ugalde, C. Tout, J. Gorres, and M. Wiescher, MNRAS 360, 375 

(2005). 

[15] J. Kuperus, Physica 31, 1603 (1965). 
[16] P. Endt, Nucl. Phys. A521, 1 (1990). 
[17] R. Azuma, private communication (2005). 

[18] G. Keyworth, P. Wilhjelm, J. G.C. Kyker, H. Newson, and E. Bilpuch, Phys. Rev. 176, 1302 
(1968). 

[19] C. Ugalde, R. Azuma, A. Couture, J. Gorres, M. Heil, K. Scheller, E. Stech, W. Tan, and 

M. Wiescher, Nucl. Phys. A758, 577c (2005). 
[20] C. Ugalde, Ph.D. thesis, University of Notre Dame (2005). 

[21] C. Ugalde, A. Couture, J. Gorres, E. Stech, and M. Wiescher, Rev. Mex. Fis. S52, 46 (2006). 
[22] C. Huan-sheng, S. Hao, Y. Fujia, and T. Jia-yong, Nucl. Instrum. Meth. B 85, 47 (1994). 
[23] J. Cseh, E. Koltay, Z. Mate, E. Somorjai, and L. Zolnai, Nucl. Phys. A413, 311 (1984). 
[24] J. Ziegler, The Stopping and Range of Ions in Matter, vol. 2-6 (Pergamon Press, 1985). 
[25] H. Lorenz-Wirzba, Ph.D. thesis, Universitat Miinster (1978). 
[26] A. Couture, private communication (2005). 
[27] A. Lane and R. Thomas, Rev. Mod. Phys. 30, 257 (1958). 
[28] C. Ugalde, in TUNL Progress Report (2006), vol. XLV, p. 38. 

[29] W. H. Press et al., Numerical recipes in C: the art of scientific computing (Cambridge Uni- 
versity Press, 1992). 

[30] S. Goriely, MOST: An updated HF model of nuclear reactions (Editions Frontires, Gif-sur- 

Yvette 1998, 1998). 
[31] L. van der Zwan and K. W. Geiger, Nucl. Phys. A284, 189 (1977). 

[32] W. A. Schier, G. P. Couchell, J. J. Egan, P. Harihar, S. C. Mathur, A. Mittler, and E. Sheldon, 

Nucl. Phys. A266, 16 (1976). 
[33] T. Rauscher, F. K. Thielemann, and K. L. Kratz, Phys. Rev. C 56, 1613 (1997). 



24 



